Edge magnetoplasmons in wide armchair graphene ribbons 
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We show that near an armchair edge of a wide graphene channel, and in the presence of a 
smooth step-like electrostatic lateral confining potential, the chirality, spectrum, spatial structure, 
and number of the fundamental edge magnetoplasmons (EMPs), in the v — 2 regime of the quantum 
Hall effect, depend strongly on the position of the Fermi level Ep. (i) When Ep is small enough and 
intersects four degenerate states of the zero Landau level (LL) at one location and two degenerate 
states of this level at a different one, two fundamental, counter propagating EMPs exist with opposite 
chirality. This is in contrast with EMPs in conventional two-dimensional electron systems in which 
only one fundamental EMP exists. For the same wave vector these EMPs have different moduli of 
phase velocities and an essential spatial overlap. These EMPs can be on resonance in a wide range 
of frequencies, for micron or submicron lengths along the edge, (ii) When Ep is sufficiently high 
and intersects only two degenerate states of the zero LL only one fundamental EMP exists with the 
usual chirality. 

PACS numbers: 71.10.Pm, 73.21.-b, 81.05.Uw 



I. INTRODUCTION 

Since the experimental discovery of graphene and the 
manifestation of its high-quality free-standing samples 1 , 
graphene has attracted a strong attention, see, e.g., Ref. 
2 for recent review. It is currently the subject of many 
independent studies because it's electronic properties are 
drastically different from those, say, of conventional two- 
dimensional electron systems (2DES) in semiconductor 
heterostructures, on liquid helium, etc. Charge carriers 
in a single-layer graphene behave like " relativistic" , chi- 
ral massless particles with a "light speed" equal to the 
Fermi velocity, vp, and possess a gapless, linear spec- 
trum close to the K and K' points-^— . One major con- 
sequence is the perfect transmission through arbitrarily 
high and wide barriers upon normal incidence, referred to 
as Klein tunneling^^, and the direction-dependent tun- 
neling through barriers^. Other unusual properties are 
a half-integer quantum Hall effect (QHE)^^— , a mini- 
mum metallic conductivity, and a zitterbewegun g 2 ' 9 ' 10 . 
The latter effect is due to a lateral confinement of Dirac 
fermions and its manifestation can be essentially modified 
by a strong magnetic field 9,10 . In addition, the submicron 
long mean-free pathsi may have important consequences 
for applications in graphene-based devices, such as tran- 
sistors, which have already been produced^. 

In addition to the aforementioned studies, graphene's 
edges have also been studied considerably- 2 -^— , in 
particular, in connection with the QHE2£~— . For some 
properties it matters whether the ed ges are of the arm- 
chair or zigzag type2£~— , see Refs. Il2lll3l for nanorib- 
bons. Also magnetic interface states (that can be un- 
derstood as a new type of edge states made within a 
" bulk" of the graphene flake) in graphene-based quantum 
wires, created by a smooth electrostatic confining poten- 
tial, have been studied 14 . There have also been many 
studies on plasmon o 15 ' 16 and magnetoplasmons 1 - , collec- 
tive wave excitations, in unconfined graphene; for a re- 



view on plasmons and magnetoplasmons in conventional 
2DES, e.g., seeA^- 19 . Different types of edge magneto- 
plasmons (EMPs) have been studied theoreticall y 18 ' 20 " — 
and experimentally^r-^ . only for conventional 2DESs, 
see, e.g., Refsi 18 ' 19 , and in particular Refs<22r— for ex- 
perimental studies of EMPs in the QHE regime. To our 
knowledge only the results of the EMP studies of Ref. I24I 
in the integral QHE regime, explain well the experimental 
findings o^ 8 - and22. for different realistic conditions. In- 
deed, the classical model of Ref. [l8| has major drawbacks 
in explaining the experimental results in the integer QHE 
regime, according to the analysis carried out in Ref. |28|. 
This is also the case with the classical model of Ref. [2l| 
according to Refs. [H, and [13. 

In this work we explore theoretically the possibility of 
fundamental EMPs in graphene following Refs. [24I25I . 
As will be shown, in the presence of a smooth yet step- 
like lateral confining potential near an armchair graphene 
edge, at y — L y /2, EMPs are possible in the v = 2 QHE 
regime. These EMPs depend strongly on the position 
of the Fermi level Ep. For case (i), referred to in the 
abstract, the main resonance (Eq. (142j) . of two EMPs of 
opposite chirality, localized within a submicron length of 
an armchair edge) is possible, e.g., if a strong coupling of 
the EMPs holds at the ends of the segment L™ % < L x , 
where L x is the length of the graphene channel. Our 
study shows that the relevant condition L x m < 1/zm is 
realistic. The EMPs that we find near a graphene edge 
are very different from those EMPs treated previously 
in conventional 2DES, in particular in the integral QHE 
regim o 18 ' 21 " — . 

In Sec. II A we present the wave functions and the 
spectra of the Landau levels (LLs) in an infinitely large 
graphene flake in the presence of a perpendicular mag- 
netic field and of a smooth electrostatic confining poten- 
tial, along the y direction. In Sec. II B we study the com- 
bined effect of a smooth, step-like electrostatic confining 
potential and of armchair graphene edges, at y = ±L y /2, 
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on the local Hall conductivity in the v = 2 QHE regime. 
In Sec. Ill we present the resulting EMPs, at an edge 
region of a wide channel, and their strong dependence on 
the position of the Fermi level Ep. We make concluding 
remarks in Sec. IV. 



II. GRAPHENE CHANNEL AND LOCAL HALL 
CONDUCTIVITY 

A. Effect of a smooth potential on the LLs 

We consider an infinitely large flat graphene flake in 
the presence of a perpendicular magnetic field B = Bz 
and of a smooth confining potential V y — V(y) along the 
y direction, of electrostatic origin. For definitcness we 
assume that this potential is symmetric. First we con- 
sider solutions with energy and wave vector close to the K 
point. In the nearest-neighbor, tight-binding model the 
one-electron Dirac Hamiltonian, for massless electrons, is 
Ti. = vpa -p-\-W v , with 1 the 2x2 unit matrix. Explicitly 
H is given by (e > 0) 

n = vp ( Jy' VF R -fit ~ eBy ) , (i) 

\p x + iPy-eBy Vy/VF 

where p x and p y are the components of the momen- 
tum operator p and vf ~ 10 6 m/s the Fermi veloc- 
ity. The vector potential is taken in the Landau gauge, 
A = (-By, 0,0). The equation {H - E)tp = admits 
solutions of the form 



^(r) = e^ x $(y)/y/L x , 



B* B (y) 



(2) 



where L x is the length of the structure along the x 
axis; due to neglected spin-dependent contributions to 
the Hamiltonian Eq. ([I} each eigenstate is two times 
degenerated on the spin quantum number. The compo- 
nents &A(y) and $3(2/) correspond to the two sublattices 
and are assumed normalized. Then the coefficients A and 
B satisfy the relation \A\ 2 + \B\ 2 = 1. To simplify the 
notation in what follows we will write &A(y) = &a and 
<1> B (y) = <1> B . Using Eqs. flTJ and © we obtain 



V y A<Z> A + hv F (k x 



where £q — (h/eB) 1 / 2 is the magnetic length. For E =^ V y 
we solve Eq. (3) for A&a and substitute the result in Eq. 
(4) . Assuming B ^ and E — V y ^ this gives 



8 2 f 2 ft 2 ?; 2 



d(E-U(£;y ))/d£ „ d 

E-U&yo) =°" (5) 



where we introduced the dimensionless variable £ = 
(y - yo)/to, w hh 7/o = Zlk x , and the notation [/(£; y ) = 
V(£ot + uo) = V(y). 

Assuming V y is a smooth function of y, with a charac- 
teristic scale Ay 3> £q, we can make the approximation 



E-U(S;yo)nE-V{yo) 



? g| k=0, (0) 



where it is used that U(0;yo) = V(yo)- 

First, for a given quantum number yo, we assume E = 
V(yo). Then we rewrite Eq. (JSJ), for £ ^ as 



a 



1 d 



(7) 



where a = dU(£; yo) / 'd£\%=o and r — £oa/hvp. Two inde- 
pendent solutions of Eq. (JT)) are given by £H*jj 2 (b£ 2 ) and 

SH^b?), with b = (i/2)(l - r 2 ) 1 / 2 , where H { yf{z) 
are the Hankel functions. The general solution is a lin- 
ear combination of them. In particular, for a constant 
electric field of arbitrary strength, we can obtain exact 
results for the wave function and eigenvalue of the n = 
LL for both r < 1 and r > 1, see Refs. HHHl 

For the sake of comparison with the well-known exact 
results of Refs. I3lll32l . for a constant electric field and 
r < 1, we first assume r < 1. Then we readily obtain 
$ fl (0 oc £Jf$(0 cx exp[-(l - r 2 ) 1/2 £ 2 /2], and 

A<J> A (£) = (l/r)[l + (l/0d/de]3$ B (O; (8) 

it follows that &a(0 = &b(0- So far the calculations 
were performed for £ ^ 0. Due to £- t/(0; yo) = 0, for 
£ = the transition from Eqs. (3) and (4) to Eqs.© 
and ([7]) is not yet justified. Here this special point can 
be dealt with easily because the continuity of the wave 
function (2J implies that of 3> A (£) and 3>b(£) at £ = 0. 

The procedure given so far applies to the K valley. If 
we repeat it for the K' valley, we obtain again Eq. (J7J. 
If we label the two valleys by k = ±, we can write both 
results in the form (— oo < £ < oo) 



^(£) = <i>B K (£) = [(l-r 2 )A£ 2 ] 



2U/4 -(l-r 2 ) 1 



A K = (l/r)[l — k(1 — 



r 2\l/2- 



2e/ \ (9) 
(10) 



y /H-8/6y)B* B = EA* A , (3) ^ ^ 



k(1 



1 

2 -j 1/2 -j -1/2 



ftwF(*x-»/^ + a/ay)A* A + ^,B*B=£7B*B > (4) and A+ = B_ 



V, 



k(1 - r 2 ) 1 / 2 ]!^ and = r(1 _ 

For r < 1 we have B+ = A- w 1 
r/2. Notice that for the linear potential 
eEy the exact results of Refs. [3l| and [32| for the 
n — LL coincide with ours, for the eigenvalue Eq ^ = 
eEy , and the wave function defined by Eqs. (9), (10). 
In what follows, we assume that r<l. 

Now, to study the effect of a smooth potential on the 
n 7^ LLs, we will use that for these LLs \E — V(j/o)| ^> 
|<z|. Then combining Eq. © with Eq. ([5]) we obtain 



— - e + J [(E - V(y )f + ] $ Bfc (£) = 0. (11) 



This is a harmonic oscillator equation whose solution is 
standard. For N = 1, 2, ... the eigenvalues are E± N ^ k;L — 
±(hv F /e )(2N) 1 / 2 + V(y ); the eigenfunctions are ap- 
proximately the well-known ones for V(y) = 0. We em- 
phasize that as here N ^ 0, the condition \E — V(yo) \ 3> 
\a\ reduces to r -C 1. 

Finally, for any n = 0, ±1,±2, ... LL and yo not too 
close to the graphene lattice termination at y = ±L y /2 
(see Fig. 1 which agrees with Ref. [lj), the eigenvalues 
iVo can be written as 



E n ,k x = sgn(n)(hv F /£ ) 



-V(y ), 



ra = 0,±l,±2,... 

(12) 

where the sign function sgn(n) = 1 and —1 for n > and 
n < 0, respectively. Notice that each n ^ LL is twice 
degenerate with respect to the valley quantum number 
k. Accordingly, for any n ^ LL and yo not too close 
to the graphene lattice termination at y = ±L y /2 (see 
Figs. 1(a), 1(b), compare with Ref. Il4h the eigenvalues 
([T2"T) are four times degenerate due to the spin and valley 
degeneracies. 



B. Effect of a smooth potential and of an armchair 
edge on LLs 

and local Hall conductivity in the v — 2 QHE regime 

Extending magnetotransport formulas for the local 
Hall conductivity cr yx (y) of a standard 2DES in the 
channel, in the presence of a smooth, lateral confining 
potential^—, we obtain, for linear responses and in 
strong magnetic fields, <J yx (y) in the forrr*^ 

<r V x(y) = n(y)e/B, (13) 
where n(y) is the y-dependent electron density given by 

n(y) = V/« re <HM(r-r)M, (14) 



with a = {n,k x }. For a finite hole density, p(y), in Eq. 
([IB"]) it follows that n(y) is changed on [n(y) —p(y)]. The 
standard relation a yx (y) = —<J xy (y) holds. Equation (fT5|) 
can be rewritten as 



°vx(y) 



4 E, 

n>0,K=± 



dyo[fn,y ,f. 



X 



\b:\ 2 \<i> 



finite, 
2 



Bk 



(Oil, (15) 



with f n ,y ,K the Fermi function and $^ K (£)> 
& Bk (C) given by Eqs. (10)-(11) in the linear-response 
limit r —> 0; the factor 2 accounts for spin degeneracy. 
Point out that k = ± in Eqs. (fl4 |) . (|15[) is actually un- 
derstood as the pseudospin quantum number; because, 
at yo > 0, only for (L y /2 — yo)/@o ^ 1 it can be well 
approximated by the valley index. Indeed, a strong split- 
ting between the electron, k = +, and the hole, k = — , 
branches of the n = LL^&ii, due to hybridization of 
the valley states take place nearby the armchair edge, 
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Figure 1: (Color online) (a) Energy spectrum of the n = 
— 1, 0, and 1 LLs as a function of the quantum number yo, 
pertinent to the right half of a symmetric graphene channel 
with armchair edges and the smooth electrostatic potential 
of finite strength, Eq. (| 16f> . for two different values of the 
Fermi level: (i) = V /2 and (ii) E ( p ] = Wo/2. For 

both cases the v = 2 quantum Hall effect is manifested in 
dc magnetotransport; here Vo = Hvf /V2£o, Ay = 10^o, and 
lo is the magnetic length. In case (i) two spatially separated 
edge states, of opposite chirality and different degeneracy, are 
created at = 5Ay as the four times degenerate n — LL 
crosses the Fermi level and at yf w 9Ay as the doubly 

degenerate branch of the n = LL crosses Ep. Both y™ 
and yf are marked by an upward arrow. In case (ii) only the 
n = LL at y"' 1 w 9Ay (downward arrow) crosses the Fermi 

(2) 

level E F , by its doubly degenerate branch that goes up with 
increasing yo. (b) The same LLs for a graphene channel in 
the absence of a smooth electrostatic potential; this spectrum 

(i) 



agrees with that of Refs. 94811131 Here for both positions E F 
and E F of the Fermi level the picture of the edge states i 
qualitatively the same as case (ii) of Fig. la. 



at | L y /2 — yo | < lo- The eigenvalues of the n = LL 
for k = +(— ) increase (decrease) with increasing yo', 
in Eq. (TT5l) it is used that the contribution from the 
(n = 0, k = — ) LL is described better by the hole rep- 
resentation. However, for the n > 1 LLs the k = ± 
branches at the armchair edge have a small splitting, due 
to hybridization of the valley states, as their eigenvalues 
increase with increasing yo', these branches are attributed 
to the electron band. 

We now consider the situations depicted in Fig. 1(a) 
for a wide symmetric channel L y > 2y r 3> Ay 3> A), 
where L y /2 — 9 Ay, y r /Ay — 5, and Ay/£ — 10. How- 
ever, y r /Ay, and Ay/lo can take any large value if the 
EMPs at the right part of the channel are well decoupled 
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from those at its left part. For clarity the smooth lateral 
potential is taken as 

V(y) = (V /2)[2+<f>((y-y r )/Ay)+$((y+y r )/Ay)\, (16) 

where §(x) is the probability integral. When the Fermi 
level Ep is between the bottoms of the n — and n = 1 
LLs, at yo = 0, and the condition Vq 3> 2kgT holds, the 
occupation of the n > 1 LLs is negligible; the same holds 
for the n = LL in the regions of yo that are well above 
Ep, see Fig. 1(a). In addition to the smoothness of the 
potential (fT6|) . we assume armchair edges of the graphene 
sheet at y = ±L y /2, which cause the bending of the 
LLs^l>i^M. In Fig. 1(a) we have L y /2-y r = 4Ay but 
our main results hold qualitatively for L y /2 — y r > Ay 
as well. 

We point out that the fine structure of the LLs n = ±1 
in Figs. 1(a), (b) , due to the removal of the pseudospin 
degeneracy at \L y /2— y Q \ < £ resulting from two possible 
hybridizations of the valley states for n ^ LLs-iI, is dis- 
carded. However, a strong splitting of the n — LL 2 ^^ 



at \L y /2 - I 



< 



due to only one possible hybridiza- 



tion of the valley states^, is taken into account. Due 
to these assumptions, for yo very close to the armchair 
termination we can formally assume that the eigenvalues 
in Eq. ([T5]) and the wave functions of the n^O LLs are 
independent of n while the eigenvalues of the n = LL, 
for k = +(— ) increase (decrease) with increasing yo- 



For either Fermi level position, E^? or Ep* 1 in Fig. 
1(a), the v — 2 quantum Hall regime will be manifested 
in dc transport measurements. Indeed, for the x axis 
normal to the plane of Fig. 1(a) and magnetic fields 
B > 0, it follows that at yo = y™ ( y" = y r = 5Ay) the 
fourfold degenerate n = LL (composed both of a dou- 
bly degenerated conduction band branch, k = +, and a 
doubly degenerated valence band branch, k = — ) crosses 
the Fermi level Ep as it goes up, with increasing yo- 
This creates a fourfold degenerate edge state that prop- 
agate along the positive x axis. However, at yo = yf 
(here yf ss 9Ay is very close to the armchair termina- 
tion of the graphene channel) only a doubly degenerated 
valence band branch (k = — , with different spin quan- 
tum numbers but the same hybridization of the valley 
states) of the n = LL goes down and crosses Ep with 
increasing y . This branch creates a doubly degenerate 
edge state that propagate along the negative x axis. We 
call this situation case (i). For the Fermi level position 



E 



(2) 

p' in Fig. 1(a) only a doubly degenerated conduction 
band branch («; = +, with different spin quantum num- 
bers and the same hybridization of the valley states) of 

(2) 

the n = LL goes up with increasing yo and crosses E F 



at yo — y" ' , where y^ 



9 Ay is very close to the arm- 



chair termination of the graphene channel. This branch 
creates a doubly degenerate edge state that propagate 
along the positive x axis. We call this situation case (ii). 

In Fig. 1(b) we plot the same LLs of the graphene 
channel in the absence of a smooth electrostatic poten- 
tial. The spectrum shown is in agreement with that of 



Refs. @-[lGJ. Here for both E ( p ] and E ( p ] the picture of 
the edge states is qualitatively the same as that for case 
(ii) of Fig. 1(a) and magnetotransport measurements will 
manifest the v = 2 QHE. Below we show that in cases 
(i) and (ii) the properties of EMPs arc very different. 

In case (i), for y > and {yf — yo)/^o ^> 1, from Eqs. 
(9)-(10), (III]), (HID and HID we obtain 



T yx {y) = — [l + exp([V%) 



4e 2 r l- 1 2e 2 

V(y?)]/k B T)] 

(17) 

where it is assumed that Eo, yo — V(yo) is smooth on 
the scale of £o, i-e., £odV(y")/dy ksT; the factor 
4 accounts for spin and pseudospin degeneracy. This 
condition of smoothness can be rewritten, upon intro- 
ducing the characteristic length It = £o(kBT£o/hVg), as 
£o <C £ti where — £$ h~ 1 dV(y^:)/dy is the group veloc- 
ity at the edge ?/". Notice that by using Vo = hvp/V2£o, 
£o/Ay = 0.1, and all other conditions applying to Fig. 1, 
we obtain = (£o/ V2tt Ay) x«fs;4x 10 6 cm/s. Also, 
for qualitatively similar conditions we obtain v^jvp = 

(£ /V2nAy) «S 1 due to the condition to/ Ay < 1. 
For sufficiently smooth potentials we can write 

W + (»"#)) * V(Vr) + (y-yr)dV(y)/dy\ y=y u, (18) 

where the second term can be written as 



{h/£l) (y - 



dE 
hdk x 



\ y ^ = (h/t 2 )vv(y-y?). (19) 



For 1 2/ — 2/t- I < £t the approximation ([TBI and Eq. (|17[) . 
for (£t/ Ay) 2 ^ 1, allow us to rewrite Eq. (fT7| as 

4e 2 r i - 1 2e 2 

*v*(v) = —[ 1 + e M(y - tf)/4r]] - —■ (20) 

We remark that setting y — y — gives 2 - 
da yx (y) 4e 2 1 



dy 



h A£ T cosh 2 (y/2£ T )' 



(21) 



Hence, in case (i), for y > and {yf — y)/£o ^ 1, we 
have Eqs. ([20]) -([2l]). Further, for L y /2 >y> L y /2 - 5£ 
we model the v = 2 numerical results ^ 13 ' 14 with the 
density 



n(y) -p(y) 



7T 3 / 2 / 3 , 

" l o J—oo 



dyoe '(y~yo) 2 /el [f 0mi _-i], 

where we assumed that Eo, yo .- is a sharply decreasing 
function at yo ~ yf such that the Fermi function in Eq. 
(|2"2"]l is very fastly growing at yo ~ yf on a scale smaller 
than £ . Then from Eqs. ([13]) and §h§ we obtain 

d<j yx (y)/dy = (2e 2 /hV^£o) e^^ 2 ^, (23) 

by changing the derivatives with respect y to those with 
respect yo and integrating by parts. 
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In a similar manner, for case (ii) and y > 0, we obtain 
that E 0tV0 + is a sharply increasing function at yo s» y"' 1 
and 

da yx (y)/dy = -(2e 2 /^o) e^"^' 1 ) 2 ^, (24) 
in agreement with Fig. 1(a). 



III. STRONG DEPENDENCE OF EMPS ON 
THE FERMI-LEVEL POSITION FOR v = 2 

Now we will study EMPs for cases (i) and (ii), see 
Fig. 1, neglecting dissipation. We expect that the charge 
excitation due to EMPs at the right part of channel will 
be strongly localized at y" (p ru (t,r)) and y d (p rd (t,r)), 
in case (i), and at y"' 1 (p r ' ul (t, r)) in case (ii). Then for 
case (i) the components of the current density j(y), in 
the low- frequency limit ui <C vf/£q, ar o 24 ' 25 

3x{y) = -<7yxE y (y) +vy ru (uj 7 k X7 y) + v d p rd (uj, k x ,y) 7 



3y(y) = Vyx{y)E x (y), 



(25) 
(26) 



where we suppressed the factor exp[— i(uit — k x x)\ com- 
mon to all terms in Eqs. (I2"5j) and ([2o| . From Eqs. (|2"5)) 
and (|26p. Poisson's equation, and the linearized continu- 
ity equation we find the integral equation for the charge 
density p(uj,k x ,y) = p ru (lu , k x , y) + p rd (uj,k x ,y) 

(u - k x v u Q )p ru {uj, k x ,y) + {uj- k x v d )p rd {uj, k x , y) 



2k x da yx (y) 
e dy 



d v' R g (\y - y'\,k x ;d) 



x [p™(u,k x ,y') + p rd (u>,k x ,y')] =0. 



(27) 



For a metallic gate placed on top of the sample, at a 
distance d from the 2DES (usually this is a heavily doped 
Si separated from the graphene sheet by a SiC>2 layer of 
thickness d = 300 nm), R g (...) is given by 

Rg{\y - y'\,k x ;d) = K (\k x \\y - y'\) 



- K (\k x W{y-y>y+Ad?l (28) 

where Kq(x) is the modified Bessel function. In the ab- 
sence of a metallic gate, d — ¥ oo, the dielectric constant 
e is spatially homogeneous if not stated otherwise. 
As da yx (y)/dy is too small according to Eqs. (f2"Tj) and 
except at y « y" and yf , we rewrite Eq. ([27)1 as 

(w - k x v%)p ru (u>, k x , y) + (u- k x v d )p rd {uj, k x ,y) 



c hk x 



_2£ T cosh 2 {y/2t T ) v^o 
dy'R g (\y - y'\,k x ;d) 

x [p^Wj + p^fofc*,!/')] = 



where = 4e 2 /he. In the long-wavelength limit 
(fe^T < 1 we have Xo(IMj/ - y')\) ~ \n{2/\k x {y - 
y')\) ~ 7; where 7 is the Euler constant. The effect of 
the gate becomes essential if d is not too large, i.e., for 
2|fc x |d < 1. For the gated sample and 4d 2 3> $to> m 
the long- wavelength limit 2|fc a |cZ -C 1, we have R g s» 
ln(2d/|y-y'|). 

From Eq. ([29]) it follows that p ru (w , k x , y) and 
p rd (u>, k x , y) can be well approximated by 



P ru {u,k x ,y) 



U T cosh z ( 



2,y-Vr 



2£ T 



) (30) 



/V fc Bj y) = (l/^o)e- fe -^ )2/£ « P r \uJ, k x ). (31) 

If we assume y d — y^^>£r, we can neglect any overlap 
between p ru (uj, k x , y) and p rd (uj, k x , y) in Eq. ([2"9"j) . Then, 
by integration of Eq. (|29l) over y within separate regions 
around y" and yj?, we obtain two coupled equations for 
p ru {uj, k x ) and p rd (uj, k x ). They read 



(w - fc^j) - 2c h k x a p (k x ; d) p ru (uj, k x ) 
2c h k x R g {\y d r ~y^\,k x ;d)p' rd {u:,k x )=Q, (32) 



with v u g = (£ /V2^Ay)v F rj 4 x 10 6 cm/s, and 

(w-MjJ + c^Om^jd) p rd (uj,k x ) 
+ c h k x R g (\yf - y%k x ;d)p ru (u,k x ) = 0, (33) 

with 3> and < 0. Indeed, we estimate a typical 
1^1 ~ 3 x 10 7 cm/s using numerical results from, e.g., 
Ref. [lj. Notice that from Refs. SHU and Sec. II we 
obtain \v d \ < vf ~ 10 8 cm/s; that is, the group velocity 
of any edge state must be smaller than Vp. The matrix 
elements a p (k x ;d) and a m (k x ;d) are given by 



a p {k x ;d) = — 



—00 «/ —00 



cosh 2 (x/2) cosh 2 (</2) 



(34) 



2 ^oo /*oo dxdt 
a m (k x ;d) = - / / _„2 +t 2 #g(41ff - 1\, k x ; d). (35) 



— OO J — OO 



(29) 



For |fc x (y^ — y")| ^> 1 it's a good approximation to 
neglect the terms oc R g (---) in Eqs. (|32|) and (|33)l . Then 
Eqs. (|32|) and (|33|) are decoupled. The resulting disper- 
sion relations for the two fundamental EMP modes are 

u+] (.kx,d) = k x Vg + 2c h k x a p {k x ;d), (36) 

for the mode localized at y", that has positive phase and 
group velocities, and 

u_ fi {k x ,d) = -k x \v d \ - c h k x a m {k x ;d), (37) 

for the mode localized at y d that has negative phase and 
group velocities. Notice that in the long- wavelength limit 
and for large d the effect of the gate, oc exp(— 2|fca;|d) <C 1, 
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k (10 3 crrT 1 ) 



100 



Figure 2: (Color online) The dispersion relations LJ^(k x , d = 

300 nm) (solid curves), u>^ (k x ,d = 300 nm) (dashed curves), 

and u£ (k x ,d — > oo) (dot-dashed curves) of two counter 
propagating fundamental EMPs for case (i) at v = 2. Panels 
(a), (b), and (c) correspond to three characteristic k x regions: 
10 3 cm _1 > k x > 10 2 cm" 1 in (a), 10 4 cm" 1 > k x > 10 3 
cm" 1 in (b), and 10 5 cm" 1 > k x > 10 4 cm" 1 in (c). The 
solid curves marked by 1, 2, and 3 correspond, respectively, 
to inter-edge distances, in the n — LL, yf — = Ay, 4Aj/, 
and 16Ay. The other parameters are B = 9T, T = 77K, 
It I la = 2, Ay = W£ , < = 4 X 10 6 cm/s, u 9 d = -3x 10 7 
cm/s, e — 2, and £o ~ 8.5 nm. 



can be neglected in Eqs. (|36p an d (|57)) . The result 
is a p (k x ;d) ->• [ln(l/|fe x |#r) - 0.145] and a m (k x ;d) -!• 
[111(1/1^14) + 3/4]. 

For case (ii), in the low- frequency limit uj <C vp/lg, 
the result is 



1 -iv-v"' 1 ) 2 /^ n r,m 



p r ^{u>,k x ). (38) 



Using Eq. (IM|) and other relevant expressions gives, for 
//'" (a;, k x ) ^ 0, the dispersion relation for only one fun- 
damental EMP mode, localized mainly at y"' 1 , 



uj {ll) {k Xl d) = k x vf +c h k x a m (k x ;d), 



(39) 



with positive phase and group velocities. Here Vg 1 > 
and, similar to \v d \, we estimate Vg 1 < 3 x 10 7 cm/s. 

If we take into account the coupling in Eqs. (I52"j) - (13111) . 
due to R g (\y?-y?\,k x ;d) ^ 0, then a nontrivial solution 
of this system requires its determinant to vanish. This 
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Figure 3: (Color online) The dispersion relations of the 
unique, left fundamental EMP for case (ii) and two values of 
d, uj {m \k x ,d = 300 nm) (solid curves), and uj (m \k x ,d = oo) 
(dashed curves). Panels (a), (b), and (c) correspond, respec- 
tively, to panels (a), (b), and (c) of Fig. 2. Here Vg 1 = 3 x 10 7 
cm/s and the other parameters are the same as those in Fig. 
2. 



leads to two renormalized EMP modes, w\/ and ui_ , 



i 



±- 



[u+fliK, d) - Lo^] (k x ,d)] 2 

1/2 



-8c 2 h k 2 x R 2 g(y d r -y^k x ;d) 



(40) 



where u)± Q (k x , d) are given by Eqs. (131))) - (13"T1) . If we 
neglect the Coulomb coupling R g (-.-) between the charge 
excitations at y d and y", Eq. (|4"0")) leads to the limits 

aj ( +\k x ,d) -t oj^] Q (k x ,d) and u>^\k x ,d) -> u)® {k x ,d). 
From Eqs. (gUJ) and ([3"2" ]) -(|3"3" ]l it follows that 

p ru (u J { l\k x ,d),k x )/p rd (u J i l\k x ,d),k x ) 

= 2p rd (u J { l\k x ,d),k x )/p ru (u J { ! ) (k x ,d),k x ), (41) 

for any d, y d — y", and k x , in particular for d — > 
oo. That is, the ratio of the charge amplitudes 
p ru {u i + ) (k x ,d),k x )/p rd (uj ( l > {k x ,d),k x ) = p r + u /pl d for the 
w+ (k x , d) EMP, at the edges y™ and y d , times the ratio 

p™lp rd , for the u { l\k x ,d) EMP, is equal to 2. 

For case (i) and v = 2, in Fig. 2 we plot the dis- 
persion relations w±(k x ,d = 300 nm) (solid curves, 
Eq. gO])), ujf> (k x ,d = 300 nm) (dashed curves, Eqs. 
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(l36|) - (j37|) L and ui±) (k x ,d — > oo) (dot-dashed curves) for 
v% = 4 x 10 6 cm/s, v d g = -3 x 10 7 cm/s, and e = 2, 
in three characteristic k x regions: 10 3 cm~ 1 > k x > 10 2 
cm -1 in (a), 10 4 cm > k x > 10 3 cm" 1 in (b), and 
10 5 cm _1 > k x > 10 4 cm -1 in (c). Here we assume that 
on one side of the graphene sheet there is SiC>2 substrate, 
with dielectric constant ~ 3, and on the other side there 
is air or vacuum: then for e we must use, in all formulas, 
an effective dielectric constant ~ 2. The other parame- 
ters used are B = 9T, T = 77K, Ay = Wt , which gives 
tr/to = 2, and to ~ 8.5 nm. The solid curves marked by 

1, 2, and 3 correspond to the inter-edge distance of the 
n = LL yf — yf = Ay, 4Ay, and 16Ay, respectively. 
For any of these curves we assume (yf — yf)/yf ^ 1; 
the dashed and dot-dashed curves are independent of 
the inter-edge distance. This allows us to neglect the 
coupling of the fundamental EMPs, localized in some re- 
gions of y > 0, with any EMPs on the left part of channel. 
In particular, for yf — yf = 16Ay the channel width is 
much larger than that used in Fig. 1. The solid curves 1, 

2, and 3 are very close in (a) to pertinent dashed curve 
and are even closer in (b) and (c). Resonances due to 
these two counter propagating EMPs (localized in a re- 
gion of extent < 1/zm, from the right edge at y yf) are 
possible in (a),(b), and (c) for L x ~ 10 _2 cm, ~ 10~ 3 cm, 
and ~ 1CP 4 cm, respectively. In Fig. 2(a) the solid and 
dashed curves show a strong effect of the gate, compare 
with the dot-dashed curves, and their behavior is very 
close to a linear one. In Fig. 2(c) both the effect of the 
gate and that of yf — yf become essentially smaller for 
the curves 1, 2, and 3. 
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Figure 4: (Color online) The ratio -p r +/p+ d = -2p r _ d /p r _ u 
versus k x for case (i), v = 2, and other conditions as in Fig. 2. 
The solid (dashed) curves correspond to d = 300 nm (d — > oo). 
The curves 1, 2, and 3 correspond to yf — y^ = Ay, 4Aj/, and 
16 Ay, respectively, and Ay = 10^o- 

For case (ii) and v = 2 in Fig. 3 we plot the dispersion 
relations of the unique left fundamental EMP (k x ,d), 
Eq. (|3T))) , for two characteristic values of d, d = 300 nm 
(solid curves) and d = oo (dashed curves), and the same 
three characteristic k x regions of Fig. 2. 

In Fig. 4, for case (i) and other conditions as in Fig. 
2, we plot the ratio —p r +lp r f versus k x . The curves 




L x em (10" 3 cm) 

Figure 5: (Color online) The main resonance frequency oj as 
function of L% m , calculated from Eq. (|42|) with N — 1, for case 
(i), v = 2, and other conditions as in Fig. 2. The resonance is 
due to two counter propagating fundamental EMPs, Eq. (|40|) . 
localized between the edge states, at y" and yf, of the n = 
LL. The solid (dashed) curves correspond to d — 300 nm 
(d — > oo). The curves 1,2, and 3 correspond to yf — y" = Ay, 
4Ay, and 16Ay, respectively, Ay = 10£o, and to = 8.5 nm. 

1, 2, and 3 are obtained from Eqs. (32J), (|40|) -(|4T1). 
The solid curves correspond to d = 300 nm and the 
dashed ones to d — > oo. According to Eq. (|41[) we have 
P+ L I ' P T + — 2p1 !i /p" 1 . Figure 4 shows that the fundamen- 
tal EMPs with positive (u>^\k x , d)/k x > 0), and negative 

phase velocity (ui_ (k x ,d)/k x < 0), renormalized by the 
inter-edge Coulomb interaction, have their charge den- 
sity amplitudes, at different edges (i.e., at y" and yf), 
in opposite phase. Moreover, for the former EMP the 
charge excitation is mainly localized at the edge y" (i.e., 
the position of edge states due to the smooth confining 
potential) as —p 7 7'lp r f > 1, whereas for the latter EMP 
it is mainly localized at the edge yf as —p™/p r _^ < 1. 
In Fig. 5, for case (i) and other conditions as in Fig. 

2, we plot the main resonance, at N = 1, obtained for 
the u ( l\k x ,d) EMPs from 

[k+(u)-k-(u)]L^ = 2irN, (42) 

where L^ n is the length of a segment of the right edge 
along which the EMPs propagate freely, see Eq. (|40[) . 
Due to the counter propagation of these two EMPs, the 
relation L^ m < L x is possible. In particular, we will 
have L| m <C L x if a strong coupling between the EMPs, 
Eq. (|40p . is introduced in the relevant high-frequency 
range within two spatial regions separated by L x m <C L x , 
along the right graphene edge(s). Actually, as we see, 
e.g., in Fig. 1(a) the right edge in case (i) consists of 
two edges, at y" and yf, with pertinent edge states (due 
to two intersections of the n = LL with the Fermi 
level) . Despite that yf is very close to the right armchair 
termination of the channel, it follows that L y /2 > yf. In 
Eq. (|4"2"| the EMP's wave vector k x (oj) is obtained from 
Eq. (|3D| abbreviated asw = u;± \k x ,d). 
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In addition, for case (i) a strong Bragg coupling is pos- 
sible due to a weak superlattice along the edge, with 
period L x m , if L x /L x m >• 1. In particular, for fre- 
quencies in the THz range, see Fig. 5, L| m < 
and L x /L x m > 10 2 correspond to rather typical lengths 
L x > 10 2 /jm in experiments. 



IV. CONCLUDING REMARKS 

We studied EMPs near an armchair edge of a wide 
graphene channel, at y = L y /2, with a smooth lateral 
potential, in the v = 2 regime of QHE and when the 
Fermi level Ep is in a gap. We showed that the position 
of Ef can strongly affect the chirality, spectrum, spatial 
structure, and the number of the fundamental EMPs. 

When E^p intersects four degenerate states of the n = 
LL at y" > and two degenerate states of this level at 
yf » yf — y" > io (case (i)), two fundamental EMPs, 
with opposite chirality, counter propagate along the right 
edge of the channel. For the same wave vector the abso- 
lute values of their phase velocities are different and they 
have spatial structure along the y axis, with an essential 
overlap in the region between the edge states, at y™ and 
y~, and their vicinity. That is, the right edge consists 
of two edges at y" and yf, with pertinent edge states, 
due to the n = LL. When the Fermi level is sufficiently 

(2) 

high, E F , and intersects only two degenerate states of 
the n — LL at y"' 1 « yf (case (ii)), only one funda- 
mental EMP exists, of the "usual" chirality for edges of 
n type conventional 2DES. 

In case (i) we found that a resonance of two fun- 
damental EMPs, of opposite chirality, on the (right) 
edge of a graphene channel, is possible in a wide re- 
gion of frequencies. The main resonance described by 
Eq. (|42|) is allowed in segment lengths L x m < L x 
along the edge. The N = 1 resonance means that the 
sum of the total increases of the wave phases of the 
to+\k x ,d) and u^}(k x ,d) EMPs, during their propaga- 
tion between the ends of L x m along the positive and neg- 
ative x axis, respectively, is equal to 2tt. This partly re- 
sembles the condition for the main resonance of a usual 
EMP, see, e.g., Refi 18 i 20 , where the EMP propagating 
along the perimeter P of a conventional 2DES, typi- 
cally P > 10~ 1 cm) 18 i 20 ' 24 , acquires a phase 2ir. More- 
over, for experimentally realistic values of L x we ob- 
tained L x m <C L x . Indeed, for frequencies in the THz 
range and L x m < 1/im the experimentally realistic val- 
ues L x ~ 10 2 /im entail L x /L x m > 10 2 . Then we can 
speculate that a strong Bragg coupling is possible due to 
a weak periodic superlattice along the edge with period 
L x m . Notice that a weak superlattice potential along the 
edge, with period > 10 2 nm, has negligible effect on a 
fundamental EMP in the QHE regime in conventional 
2DESs 27 . In addition, as two renormalized fundamental 
EMP modes uf* , Eq. gO), are counter propagating and 
their spatial structures have essential overlap along y, in 



a narrow region between y" and yf, that can easily be 
modified due to the strong dependence of y" on a smooth 
lateral electrostatic potential, time-resolved experiments, 
such as those of Refs. I2II29I , can be used to observe ap- 
pearance of counter propagating EMP along the armchair 
edge in the v = 2 QHE regime. As far as, shown in Fig. 
1(a), case (i) is realized, that can be realized for a wide 
range of parameters. The entire EMP picture and prop- 
erties are different from those of EMPs in conventional 
2DESs due to the difference in the spectrum of the edges 
and the corresponding wave functions. In fact, as was 
mentioned earlier in Fig. 1(b) and above Eq. (17), in the 
absence of a smooth electrostatic potential the spectrum 
agrees with the usual, hard- wall potential of Refs. [H-HIlI. 
Correspondingly, we don't have two counterpropagating 
fundamental EMPs but only one fundamental EMP with 
properties similar to those of the fundamental EMP in 
conventional 2DES. 

Next we list and discuss the approximations used. In 
studying the EMPs in the v — 2 QHE regime we ne- 
glected dissipation. This approximation is well justified 
as the EMP damping can be related only with inelastic 
scattering processes within narrow temperature belts, of 
width fcflT, of each edge state, cfi 24 ' 25 , that are much 
weaker than scattering processes due to a static disor- 
der, especially in the QHE regime which implies relatively 
low T. The latter makes a dominant contribution to the 
transport scattering time in a 2DES of graphene 1 -! 2 -^ for 
B = 0. Further, the damping of the EMPs will influence 
some properties of a Bragg coupling and the quality of 
the EMP resonances, cf. Eq. (|%2"j) . Notice that for de- 
creasing temperature T any EMP damping will quickly 
weaken. However, for sufficiently small T the condition 
£t/£o 1 can be violated. Nevertheless, at quite low 
T and for sufficiently smooth bare confining potential, 
the group velocity can essentially decrease with decreas- 
ing temperature 2 ^ due to many-body effects. In addi- 
tion, even for £t/£o < 1 it appears that the present 
results will be only weakly and quantitatively modified 
since the contributions to a fundamental EMP coming 
from a region of the LL edge, at y", will bring about 
only small changes^—. Obviously, for a more accurate 
account of the EMPs studied here dissipation must be in- 
cluded in the treatment. We emphasize that our study of 
the fundamental EMPs for the armchair termination of a 
graphene channel cannot be directly extended to zigzag 
termination as some important properties of the wave 
functions and the energy levels are different than those 
of the armchair termination, cfj 2 -^— We relegate the 
study of EMPs along zigzag edges to a future work. 

We used a simple analytical model of a smooth, lat- 
eral confining potential Eq. (1161) . but our main results 
are quite robust to modifications of its form and param- 
eters since cases (i) and (ii) can be realized in a graphene 
channel in the v = 2 QHE regime. Further, near the 
edge states at yf and y"' 1 we used a simple analytic 
model to approximate a static density profile, cf. Eq. 
(f2"2")l . In doing so we neglected possible modifications 
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of the static density profile due to local chargin g 34 ' 38 ' 39 
oc d 2 V{y)/dy 2 . Notice that these modifications are weak 
for a smooth potential and their neglect should have a mi- 
nor effect in the fundamental EMPs wc studied as their 
main properties are very robust against details of a static 
density profile^— and, in particular, nonlocal effects^. 
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